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Alignment of the ALICE Inner Tracking System with 
cosmic-ray tracks 



ALICE collaboration 

Abstract: ALICE (A Large Ion Collider Experiment) is the LHC (Large Hadron Collider) exper- 
iment devoted to investigating the strongly interacting matter created in nucleus-nucleus collisions 
at the LHC energies. The ALICE ITS, Inner Tracking System, consists of six cylindrical layers 
of silicon detectors with three different technologies; in the outward direction: two layers of pixel 
detectors, two layers each of drift, and strip detectors. The number of parameters to be determined 
in the spatial alignment of the 2198 sensor modules of the ITS is about 13,000. The target align- 
ment precision is well below 10 /xm in some cases (pixels). The sources of alignment information 
include survey measurements, and the reconstructed tracks from cosmic rays and from proton- 
proton collisions. The main track-based alignment method uses the Millepede global approach. 
An iterative local method was developed and used as well. We present the results obtained for the 
ITS alignment using about 10^ chai^ged tracks from cosmic rays that have been collected during 
summer 2008, with the ALICE solenoidal magnet switched off. 

Keywords: Particle tracking detectors (Solid-state detectors); Detector alignment and 
calibration methods. 
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1. Introduction 

The ALICE experiment [1] will study nucleus-nucleus, proton-proton and proton-nucleus colli- 
sions at the CERN Large Hadron Collider (LHC). The main physics goal of the experiment is to 
investigate the properties of strongly-interacting matter in the conditions of high energy density 
(> 10 GeV/fm^) and high temperature ( > 0.2 GeV), expected to be reached in central Pb-Pb col- 
lisions at y/s^ = 5.5 TeV. Under these conditions, according to lattice QCD calculations, quark 
confinement into colourless hadrons should be removed and a deconfined Quark-Gluon Plasma 
should be formed [2]. In the past two decades, experiments at CERN-SPS (y^s^^ = 17.3 GeV) 
and BNL-RHIC (y^^ = 200 GeV) have gathered ample evidence for the formation of this state of 
matter [3]. 

The ALICE experimental apparatus, shown in Fig. |I], consists of a central barrel, a forward 
muon spectrometer and a set of small detectors in the forward regions for trigger and other func- 
tions. The coverage of the central barrel detectors allows the tracking of particles emitted within a 
pseudo-rapidity range |t] | < 0.9 over the full azimuth. The central ban^el is sun^ounded by the lai'ge 
L3 magnet that provides a field B = 0.5 T. 

The ITS (Inner Tracking System) is a cylindrically-shaped silicon tracker that surrounds the 
interaction region. It consists of six layers, with radii between 3.9 cm and 43.0 cm, covering the 
pseudo-rapidity range |t7| < 0.9. The two innermost layers are equipped with Silicon Pixel De- 
tectors (SPD), the two intermediate layers contain Silicon Drift Detectors (SDD), while Silicon 
Strip Detectors (SSD) are used on the two outermost layers. The main task of the ITS is to pro- 
vide precise track and vertex reconstruction close to the interaction point. In particular, the ITS 
was designed with the aim to improve the position, angle, and momentum resolution for tracks 




Figure 1. General layout of the ALICE experiment [1]. 
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reconstructed in the Time Projection Chamber (TPC), to identify the secondary vertices from the 
decay of hyperons and heavy flavoured hadrons, to reconstruct the interaction vertex with a reso- 
lution better than 100 /im, and to recover particles that are missed by the TPC due to acceptance 
limitations (very low momentum particles not reaching the TPC and very high momentum ones 
propagating along the 10% inactive area between adjacent TPC chambers). 

The measurement of charm and beauty hadron production in Pb-Pb collisions at the LHC is 
one of the main items of the ALICE physics program, because it will allow to investigate the mech- 
anisms of heavy-quark propagation and hadronization in the large, hot and dense medium formed 
in high-energy heavy-ion collisions and it will serve as a reference for the study of the effects of the 
medium on quarkonia states [4]. To measure the separation, from the interaction vertex, of the de- 
cay vertices of heavy flavoured hadrons, which have mean proper decay lengths cT ~ 100-500 /xm, 
requires a resolution on the track impact parameter (distance of closest approach to the vertex) well 
below 100 /im. This requirement is met by the ITS. The design position resolution in the plane 
transverse to the beam line for charged-pion tracks reconstructed in the TPC and in the ITS is ex- 
pected to be approximately 10 /im + 53 lJ,m/{ptVsind), where pt is the transverse momentum in 
GeV/c and 6 is the polar angle with respect to the beam line [4]. The ITS is made of thousands 
of separate modules, whose position is different from the ideal due to the limitations associated 
with the assembly and integration of the different components, and the forces these components 
experience. In order to achieve the required high precision on the track parameters, the relative 
position (location and orientation) of every module needs to be determined precisely. We refer to 
the procedure used to determine the modules relative position as alignment. The ITS alignment 
procedure starts from the positioning survey measurements performed during the assembly, and 
is refined using tracks from cosmic-ray muons and from particles produced in LHC pp collisions. 
Two independent methods, based on tracks-to-measured-points residuals minimization, are con- 
sidered. The first method uses the Millepede approach [5], where a global fit to all residuals is 
performed, extracting all the alignment parameters simultaneously. The second method performs 
a (local) minimization for each single module and accounts for correlations between modules by 
iterating the procedure until convergence is reached. 

In this article, we present the alignment methods for the ITS and the results obtained using the 
cosmic-data sample collected during summer 2008 with B = (a small data set with B = ±0.5 T 
was also collected; we used it for a few specific validation checks). In section we describe in 
detail the ITS detector layout and in section we discuss the strategy adopted for the alignment. 
In section ^ we describe the 2008 sample of cosmic-muon data. These data were used to validate 
the available survey measurements (section^ and to apply the track-based alignment algorithms: 
the Millepede method (section ^ and a local method that we are developing (section [7|). We draw 
conclusions in section S 

2. ITS detector layout 

The geometrical layout of the ITS layers is shown in the left-hand panel of Fig. 0, as it is im- 
plemented in the ALICE simulation and reconstruction software framework (AliRoot [6]). The 
ALICE global reference system has the z axis on the beam line, the x axis in the LHC (horizontal) 
plane, pointing to the centre of the accelerator, and the y axis pointing upward. The axis of the ITS 
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Figure 2. Layout of the ITS (left) and orientation of the ALICE global (middle) and ITS -module local (right) 
reference systems. The global reference system has indeed its origin in the middle of the ITS, so that the z 
direction coincides with the beam line. 

barrel coincides with the z axis. The module local reference system (Fig. 0, right) is defined with 
the x\oc and zioc axes on the sensor plane and with the zioc axis in the same direction as the global 
z axis. The local x direction is approximately equivalent to the global r(p. The alignment degrees 
of freedom of the module are translations in x\oc, Jioo Zioo and rotations by angles i/Zjoc, 0ioc> ^loo 
about the ;cioc, Jioo Zioc axes, respectively ^ 

The ITS geometry in AliRoot is described in full detail, down to the level of all mechani- 
cal structures and single electronic components, using the ROOT [7] geometrical modeler. This 
detailed geometry is used in Monte Carlo simulations and in the track reconstruction procedures, 
thereby accounting for the exact position of the sensor modules and of all the passive material that 
determine particle scattering and energy loss. 

The geometrical parameters of the layers (radial position, length along beam axis, number of 
modules, spatial resolution, and material budget) are summarized in Table |I]. The material budget 
reported in the table takes into account the -averaged material (including the sensors, electronics, 
cabling, support structures, and cooling) associated with radial paths through each layer. Another 
L30% of radiation length comes from the thermal shields and supports installed between SPD 
and SDD barrels and between SDD and SSD barrels, thus making the total material budget for 
perpendicular tracks equal to 7.66% of Xq. 

In the following, the features of each of the three sub-detectors (SPD, SDD and SSD) that are 
relevant for alignment issues are described (for more details see [1]). 

2.1 Silicon Pixel Detector (SPD) 

The basic building block of the ALICE SPD is a module consisting of a two-dimensional sensor 
matrix of reverse-biased silicon detector diodes bump-bonded to 5 front-end chips. The sensor 

'The alignment transformation can be expressed equivalently in terms of the local or global coordinates. 
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Table 1. Characteristics of the six ITS layers. 
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Material 
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modules 


rep X z [mm^] 


r(p X z LiUm ] 


X/Xq [%] 


1 


pixel 


3.9 


14.1 


80 


12.8x70.7 


12x100 


1.14 


2 


pixel 


7.6 


14.1 


160 


12.8x70.7 


12x100 


1.14 


3 


drift 


15.0 


22.2 


84 


70.17x75.26 


35x25 


1.13 


4 


drift 


23.9 


29.7 


176 


70.17x75.26 


35x25 


1.26 


5 


strip 


38.0 


43.1 


748 


73x40 


20x830 


0.83 


6 


strip 


43.0 


48.9 


950 


73x40 


20x830 


0.86 



matrix consists of 256 x 160 cells, each measuring 50 fxm (rep) by 425 /xm (z). 

Two modules are mounted together along the z direction to form a 141.6 mm long half-stave. 
Two mirrored half-staves are attached, head-to-head along the z direction, to a carbon-fibre support 
sector, which also provides cooling. Each sector (see Fig. 0, right) supports six staves: two on the 
inner layer and four on the outer layer. The sensors are mounted in such a way that there is a 2% 
overlap between the active regions in r(p, but along z there is a gap between each two consecutive 
sensors. Five sectors are then mounted together to form a half-barrel and finally the two (top and 
bottom) half-barrels are mounted around the beam pipe to close the full barrel (shown in the left- 
hand side of Fig. 0), which is actually composed of 10 sectors. In total, the SPD includes 60 staves, 
consisting of 240 modules with 1200 readout chips for a total of 9.8 x 10^ cells. 

The spatial precision of the SPD sensor is determined by the pixel cell size, the track inci- 
dence angle on the detector, and by the threshold applied in the readout electronics. The values of 
resolution along rep and z extracted from beam tests are 12 and 100 /^m, respectively. 




Figure 3. SPD drawings. Left: the SPD barrel and the beam pipe (radius in mm). Right: a Carbon Fibre 
Support Sector. 
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2.2 Silicon Drift Detector (SDD) 



The basic building block of the ALICE SDD [8] is a module divided into two drift regions where 
electrons move in opposite directions under a drift field of w 500 V/cm (see Fig. ^ right), with 
hybrids housing the front-end electronics on either side. The SDD modules are mounted on linear 
structures called ladders. There are 14 ladders with six modules each on the inner SDD layer 
(layer 3), and 22 ladders with eight modules each on the outer SDD layer (layer 4). Modules and 
ladders are assembled to have an overlap of the sensitive areas larger than 580 /im in both rep and z 
directions, so as to provide full angular coverage over the pseudo-rapidity range |t} | < 0.9 (Fig. ||, 
left). 

The modules are attached to the ladder space frame and have their anode rows parallel to the 
ladder axis (z). The ladders are mounted on a support structure made of two cones and four support 
rings to form the two cylindrical layers [9]. The support rings aie mechanically fixed to the cones 
and bear ruby spheres, used as a reference for the ladder positioning as well as for the geometrical 
survey of the module positions in the ladder reference system. 

The z coordinate is reconstructed from the centroid of the collected charge along the anodes. 
The position along the drift coordinate (xioc ~ r(p) is reconstructed starting from the measured drift 
time with respect to the trigger time. An unbiased reconstruction of the ;cioc coordinate requires 
therefore to know with good precision the drift velocity and the time-zero (to), which is the mea- 
sured drift time for particles with zero drift distance. The drift velocity depends on temperature (as 
r^^^) and it is therefore sensitive to temperature gradients in the SDD volume and to temperature 
variations with time. Hence, it is important to calibrate this parameter frequently during the data 
taking. Three rows of 33 MOS charge injectors are implanted at known distances from the collec- 
tion anodes in each of the two drift regions of a SDD module [10] for this purpose, as sketched 
in Fig. (right). Finally, a correction for non-uniformity of the drift field (due to non-linearities 
in the voltage divider and, for a few modules, also due to significant inhomogeneities in dopant 
concentration) has to be applied. This correction is extracted from measurements of the systematic 
deviations between charge injection position and reconstructed coordinates that was performed on 




Figure 4. Left: scheme of the SDD layers. Right: scheme of a SDD module, where the drift direction is 
parallel to the .:>cioc coordinate. Units are millimeters. 



all the 260 SDD modules with an infrared laser [11]. 

The spatial precision of the SDD detectors, as obtained during beam tests of full-size pro- 
totypes, is on average 35 /xm along the drift direction ^cioc and 25 /im for the anode coordinate 

Zloc- 

2.3 Silicon Strip Detector (SSD) 

The basic building block of the ALICE SSD is a module composed of one double-sided strip 
detector connected to two hybrids hosting the front-end electronics. Each sensor has 768 strips on 
each side with a pitch of 95 /^m. The stereo angle is 35 mrad, which is a compromise between 
stereo view and reduction of ambiguities resulting from high particle densities. The strips are 
almost parallel to the beam axis (z-direction), to provide the best resolution in the r(p direction. 

The modules are assembled on ladders of the same design as those supporting the SDD [9]. 
A view of the SSD ladder is shown in Fig. 0. The innermost SSD layer (layer 5) is composed of 
34 ladders, each of them being a linear array of 22 modules along the beam direction. Layer 6 
(the outermost ITS layer) consists of 38 ladders, each made of 25 modules. In order to obtain full 
pseudo-rapidity coverage, the modules ai^e mounted on the ladders with small overlaps between 
successive modules, that are staggered by 600 /im in the radial direction. The 72 ladders, carrying 
a total of 1698 modules, are mounted on support cones in two cylinders. Carbon fiber is lightweight 
(to minimize the interactions) and at the same time it is a stiff material allowing to minimize the 
bending due to gravity, which is expected to give shifts of at most 50 /im, for the modules at the 
centre of the lateral ladders of the outer SSD layer. 

For each layer, neighbouring ladders are mounted at one of two sUghtly different radii (Ar = 
6 mm) such that full azimuthal coverage is obtained. The acceptance overlaps, present both along 
z and r(p, amount to 2% of the SSD sensor surface. The positions of the sensors with respect to 
reference points on the ladder were measured during the detector construction phase, as well as the 
ones of the ladders with respect to the support cones. 

The spatial resolution of the SSD system is determined by the 95 /im pitch of the sensor 
readout strips and the charge-sharing between those strips. Without making use of the analogue 
information the r.m.s spatial resolution is 27 jUm. Beam tests [12] have shown that a spatial resolu- 
tion of better than 20 /im in the r<p direction can be obtained by analyzing the charge distribution 
within each cluster. In the direction along the beam, the spatial resolution is of about 830 /im. 

The SSD gain calibration has two components: overall calibration of ADC values to energy 
loss and relative calibration of the P and N sides. This charge matching is a strong point of double 
sided silicon sensors and helps to remove fake clusters. Both the overall and relative calibration 




Figure 5. View of one SSD ladder (from layer 5) as described in the AliRoot geometry. 
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are obtained from the data. Since the signal-to-noise ratio is larger than 20, the detection efficiency 
does not depend much on the details of the gain calibration. 

3. Alignment target and strategy 

For silicon tracking detectors, the typical target of the alignment procedures is to achieve a level 
of precision and accuracy such that the resolution on the reconstructed track parameters (in par- 
ticular, the impact parameter and the curvature, which measures the transverse momentum) is not 
degraded significantly with respect to the resolution expected in case of the ideal geometry without 
misalignment. For the ALICE ITS, this maximum acceptable degradation has been conventionally 
set to 20% (a similar target is adopted also for the ATLAS Inner Detector [13]). The resolutions 
on the track impact parameter and curvature are both proportional to the space point resolution, in 
the limit of negligible multiple scattering effect (large momentum). If the residual misalignment is 
assumed to be equivalent to random gaussian spreads in the six alignment parameters of the sensor 
modules, on which space points are measured, a 20% degradation in the effective space point res- 
olution (hence 20% degradation of the track parameters in the large momentum limit) is obtained 
when the misalignment spread in a given direction is V 120%^ — 100%^ « 70% of the intrinsic 
sensor resolution along that direction. With reference to the intrinsic precisions listed in Table |T], 
the tai^get residual misalignment spreads in the local coordinates on the sensor plane are: for SPD, 
8 /xm in xioc and 70 /im in zioc; for SDD, 25 /im in x\oc and 18 /^m in zioc', for SSD, 14 /im in 
xioc and 500 /im in zioc- Note that these spreads represent effective alignment spreads, including 
the significant effect of the dioc angle (rotation about the axis normal to the sensor plane) on the 
spatial resolution. In any case, these target numbers are only an indication of the precision that is 
required to reach an acceptable alignment quality. We will aim at getting even closer to the design 
performance expected in case of ideal geometry. 

The other alignment parameters (yioc, Yioc, (fioc) describe movements of the modules mainly 
in the radial direction. These have a small impact on the effective resolution, for tracks with a 
small angle with respect to the normal to the module plane, a typical case for tracks coming from 
the interaction region. However, they are related to the so-called weak modes: correlated misalign- 
ments of the different modules that do not affect the reconstructed tracks fit quality (x^), but bias 
systematically the track parameters. A typical example is radial expansion or compression of all 
the layers, which biases the measured track curvature, hence the momentum estimate. Con^elated 
misalignments for the parameters on the sensor plane (xioc, Zioo ^loc) can determine weak modes 
as well. These misalignments ai^e, by definition, difficult to determine with tracks from collisions, 
but can be addressed using physical observables [14] (e.g. looking for shift in invariant masses of 
reconstructed decay particles) and cosmic-ray tracks. These offer a unique possibility to correlate 
modules that are never coiTclated when using tracks from the interaction region, and they offer a 
broad range of track-to-module-plane incidence angles that help to constrain also the yioc, Wioc and 
<Pioc pai^ameters, thus improving the sensitivity to weak modes. 

As already mentioned in the introduction, the sources of alignment information that we use 
are the survey measurements and the reconstructed space points from cosmic-ray and collision 
particles. These points are the input for the software alignment methods, based on global or local 
minimization of the residuals. 
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The general strategy for the ITS first aUgnment starts with the vaUdation of the construction 
survey measurements of the SSD detector with cosmic-ray tracks and continues with the software 
alignment of the SPD and the SSD detectors, which also uses cosmic-ray tracks, collected without 
magnetic field. The initial alignment is more robust if performed with straight tracks (no field), 
which help to avoid possible biases that can be introduced when working with curved tracks (e.g. 
radial layer compression/expansion). Then, the already aligned SPD and SSD are used to confirm 
and refine the initial time-zero calibration of SDD, obtained with SDD standalone methods. These 
first steps are described in this report, which presents the status of the ITS alignment before the 
start of the LHC with proton-proton collisions. 

The next step will be, after the validation of the SDD survey measurements with cosmic-ray 
tracks, the aUgimient of the full detector (SPD, SDD, SSD) with tracks from cosmic rays and, 
mainly, from proton-proton collisions collected with magnetic field B = and B = 0.5 T. In 
particular, the data with magnetic field switched on will allow us to study the track quality and 
precision as a function of the measured track momentum, thus separating the detector resolution 
and residual misaUgnment from multiple scattering. The tracks from collisions will provide a 
uniform coverage of the detector modules and will also be used to routinely monitor the quality 
of the alignment during data taking, and refine the corrections if needed. The last step will be 
the relative alignment of the ITS and the TPC with tracks, when both detectors will be internally 
aligned and calibrated. In addition, the relative movement of the ITS with respect to the TPC is 
being monitored, and to some extent measured, using a dedicated system based on lasers, mirrors 
and cameras [15]. 

4. Cosmic-ray run 2008: data taking and reconstruction 

During the 2008 cosmic run, extending from June to October, about 10^ events with reconstructed 
tracks in the ITS were collected. In order to simpHly the first alignment round, the solenoidal 
magnetic field was switched off during most of this data taking period. The status of the three ITS 
sub-detectors during the data taking is summarized in the following paragraph (for more details 
on the sub-detectors commissioning, see Refs. [16-18]). The corresponding status during the first 
LHC runs with proton-proton collisions is given in Ref. [19]. 

For the SPD, 212 out of 240 modules (88%) were active. Noisy pixels, corresponding to less 
than 0.15% of the total number of pixels, were masked out, and the information was stored in the 
Offline Conditions Database (OCDB) to be used in the offline reconstruction. For the SDD, 246 out 
of 260 modules (95%) participated in the data acquisition. The baseline, gain and noise for each 
of the 133,000 anodes were measured every 24 hours by means of dedicated calibration runs that 
allowed us also to tag noisy (?a 0.5%) and dead (1%) channels. The drift velocities were measured 
with dedicated injector runs collected every 6 hours, stored in the OCDB and successively used in 
the reconstruction. For the SSD, 1477 out of 1698 modules (87%) were active. The fraction of bad 
strips was 1 .5%. The normaUzed difference in P- and N-charge had a FWHM of 1 1%. The gains 
proved to be stable during the data taking. 

The events to be used for the ITS alignment were collected with a trigger provided by the pixel 
detectors (SPD). The SPD FastOR trigger [1] is based on a programmable hit pattern recognition 
system (on FPGA) at the level of individual readout chips (1200 in total, each reading a sensor area 
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of about 1.4 X 1.4 cm^). This trigger system allows for a flexible selection of events of interest, 
for example high-multiplicity proton-proton collisions, foreseen to be studied in the scope of the 
ALICE physics program. For the 2008 cosmic run, the trigger logic consisted of selecting events 
with at least one hit on the upper half of the outer SPD layer (r ^ 7 cm) and at least one on 
the lower half of the same layer. This trigger condition enhances significantly the probability of 
selecting events in which a cosmic muon, coming from above (the dominant component of the 
cosmic -ray particles reaching the ALICE cavern placed below 30 m of molasse), traverses the 
full ITS detector. This FastOR trigger is very efficient (more than 99%) and has purity (fraction of 
events with a reconstructed track having points in both SPD layers) reaching about 30^0%, limited 
mainly by the radius of the inner layer 4 cm) because the trigger assures only the passage of 
a particle through the outer layer (?» 7 cm). For the FastOR trigger, typically 77% of the chips 
(i.e. about 90% of the active modules) could be configured and used. The trigger rate was about 
0.18 Hz. 

The following procedure, fully integrated in the AliRoot framework [6], is used for track 
reconstruction. After the cluster finding in the ITS (hereafter, we will refer to the clusters as 
"points"), a pseudo primary vertex is created using three aligned points in two consecutive layers 
(starting the search from the SPD). Track reconstruction is then performed using the ITS standalone 
tracker (as described in [4,20]), which finds tracks in the outward direction, from the innermost 
SPD layer to the outermost SSD layer, using the previously found pseudo primary vertex as its seed; 
all found tracks are then refitted using the standard Kalman-filter fit procedure as implemented in 
the default ITS tracker. During the track refit stage, when the already identified ITS points are used 
in the Kalman-filter fit in the inward direction, in order to obtain the track parameters estimate at 
the (pseudo) vertex, "extra" points are searched for in the ITS module overlaps. For each layer, a 
search road for these overlap points in the neighbouring modules is defined with a size of about 
seven times the current track position error. Currently, the "extra" points are not used to update the 
track parameters, so they can be exploited as a powerful tool to evaluate the ITS alignment quality. 

A clean cosmic event consists of two separate tracks, one "incoming" in the top part of the ITS 
and one "outgoing" in the bottom part. Their matching at the reference median plane (y = 0) can 
be used as another alignment quality check. These two track halves are merged together in a single 
array of track points, which is the single-event input for the track-based alignment algorithms. A 
typical event of this type, as visualized in the ALICE event display, is shown in Fig. ^ 

The uncorrected zenith-azimuth 2D distribution of the (merged) tracks with at least eight 
points in the ITS is shown in Fig. [7], where the azimuth angle is defined in a horizontal plane start- 
ing from the positive side of the z global axis. The modulations in the azimuthal dependence of the 
observed flux are due to the presence of inhomogeneities in the molasse above the ALICE cavern, 
mainly the presence of two access shafts. These are seen as the structures at zenith angle 30° 
and azimuth w 180° (large shaft) and 270° (small shaft). On top of these structures, the effect 
of the SPD outer layer geometrical acceptance is visible: the azimuthal directions perpendicular to 
the z axis (around 90° and 270°) have larger acceptance in the zenith angle. 

The main limitation of the usage of cosmic-ray tracks for the alignment of a cylindrical de- 
tector like the ITS is that the occupancy of the side modules (zenith angles approaching 90°) is 
small, especially for the external layers [21]. In the case of the SSD outer layer, which has the 
smallest fractional coverage, about 75% of the ladders are covered. This is due to the small size 
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Figure 6. (colour online) A clean cosmic event reconstructed in the ITS (left), as visualized in the ALICE 
event display. The zoom on the SPD (right) shows an "extra" point in one of the r(p acceptance overlaps of 
the outer layer 




azimuth angle [deg] 

Figure 7. (colour online) Uncorrected distribution of the zenith-azimuth angles of the cosmic tracks recon- 
structed in the ITS. 



of the triggering detector (SPD), the dominance of small zenith angles for cosmic -ray particles and 
the cut on the track-to-module incidence angle (> 30°) that we apply to reject large and elongated 
clusters. 

5. Validation of the survey measurements with cosmic-ray tracks 

The SDD and SSD were surveyed during the assembling phase using a measuring machine. The 
survey, very similar for the two detectors, was carried out in two stages: the measurement of the 
positions of the modules on the ladders and the measurement of the positions of the ladder end 
points on the support cone. 
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In the first stage, for SDD for example [22], the three-dimensional positions of six reference 
markers engraved on the detector surface were measured for each module with respect to ruby 
reference spheres fixed to the support structure. The precision of the measuring machine was 5 jxra 
in the coordinates on the ladder plane and about 10 /im in the direction orthogonal to the plane. 
The deviations of the reference marker coordinates on the plane with respect to design positions 
showed an average value of 1 jxra and a r.m.s. of 20 /im. In the second stage, the positions of 
the ladder end points with respect to the cone support structure were measured with a precision 
of about 10 nm. However, for the outer SSD layer, the supports were dismounted and remounted 
after the survey; the precision of the remounting procedure is estimated to be around 20 /im in the 
r(p direction [1]. 

In the following we describe the results for the validation of the SSD survey measurements 
with cosmic-ray data. The validation of the SDD survey will be performed after completion of the 
detector calibration. 

5.1 Double points in SSD module overlaps 

As already mentioned, the modules are mounted with a small (2 mm) overlap for both the longi- 
tudinal (z, modules on the same ladder) and transverse directions (rep, adjacent ladders). These 
overlaps allow us to verify the relative position of neighbouring modules using double points pro- 
duced by the same particle on the two modules. Since the two points are very close in space and 
the amount of material crossed by the particle between the two points is very limited, multiple 
scattering can be neglected. 

We define the distance Axioc between the two points in the local x direction on the module plane 
(f« rep) by projecting, along the track direction, the point of one of the two modules on the other 
module plane. Figure |8| (left) shows the Axioc distribution with and without the survey corrections, 
for both SSD layers (it was verified that the distributions for the two layers are compatible [21, 
23]) . When the survey corrections are applied, the spread of the distributions, obtained from a 
gaussian fit, is a f« 25.5 /im. This arises from the combined spread of the two points, thus the 
con^esponding effective position resolution for a single point is estimated to be smaller by a factor 
l/\/2, i.e. « 18 /im, which is compatible with the expected intrinsic spatial resolution of about 
20 /^m. This indicates that the residual misalignment after applying the survey is negligible with 
respect to the intrinsic spatial resolution. This validation procedure was confirmed using Monte 
Carlo simulations of cosmic muons in the detector without misalignment, which give a spread in 
Axioc of about 25 /^m, in agreement with that obtained from the data. 

5.2 Track-to-point residuals in SSD 

Another test that was performed uses two points in the outer SSD layer to define a straight track 
(no magnetic field) and inspects the residuals between points on the inner layer and the track. The 
residuals are calculated using the position along the track corresponding to the minimum of the 
weighted (dimensionless) distance to the point^. Figure |8| (right) shows the distribution of the r(p 

^The different expected resolutions in r(p and z iiave been taken into account in the calculation of the distance of 
closest approach by dividing the deviations by the expected uncertainties, i.e. making use of a dimensionless distance 
measure. 
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Figure 8. (colour online) SSD survey validation. Left: distribution of A.ji:ioc, the distance between two points 
in the module overlap regions along z on the same ladder. Right: distribution of the r<p residuals between 
straight-line tracks defined from two points on layer 6 and the corresponding points on layer 5. In both cases, 
gaussian fits to the distributions with survey applied are shown (in right-hand panel the fit range is ±2 a, i.e. 
[-60 jUm,+60 ^m]). 



residuals between tracks through layer 6 and points on layer 5. The distribution exhibits significant 
non-gaussian tails, due to multiple scattering of low-momentum particles. The effect of multiple 
scattering on the residuals was analytically estimated to be of about 300 pLra for p = 0.5 GeV/c 
and to be negligible for /? > 7 GeV/c. This is roughly compatible with the observed distribution of 
residuals. The width of the central part of the distribution is quantified by performing a gaussian fit 
truncated at 2 a, that gives Or(p — 29 IXTSX. The spread contains a contribution from the uncertainty 
in the track trajectory due to the uncertainties in the points on the outer layer. Assuming the same 
resolution on the outer and inner layers and taking into account the geometry of the detector, the 
effective single point resolution spread is l/\/l.902 times the overall spread [23], that is 21 /xm. 
This spread is larger than the effective resolution of about 18 pLra that is extracted from the double 
points in module overlaps. This difference could be partly due to the multiple scattering, relevant 
for this analysis and negligible for the overlaps analysis, but we can not rule out that additional 
misalignments with a r.m.s. up to about 15 ;Um are present in the SSD. The mean residual is also 
non-zero, (3.7 ±0.4) /xm, which suggests that residual shifts at the 5-10 pLra level could be present. 
These misalignments would have to be at the ladder level to be compatible with the result from the 
study with sensor module overlaps. 

The same analysis was performed for the residuals in the z direction [21, 23], not shown here. 
The distributions without and with survey were found to be compatible and the corresponding 
effective single point resolution was found to be compatible with the expected intrinsic resolution 
of about 800 /im. This indicates that the residual misalignment in z is much smaller than the 
intrinsic SSD resolution. 

A third method that was used to verify the SSD survey consisted in performing tracking with 



-21 - 



pairs of points (2 points on layer 5 and two points on layer 6 or two sets of points on layer 5 and 
6), and comparing the track parameters of both track segments. The conclusion from this method 
is consistent with the results from the track-to-point method. For details see [23] . 

6. ITS alignment with Millepede 

In general, the task of the track-based alignment algorithms is determining the set of geometry 
parameters that minimize the global of the track-to-point residuals: 

;tglobal= L ^t'^p^7,p5t,p. (6.1) 
modules, tracks 

In this expression, the sum runs over all the detector modules and all the tracks in a given data set; 
St,p = rt — rp is the residual between the data point Tp and the reconstructed track extrapolation r, 
to the module plane; Vf p is the covariance matrix of the residual. Note that, in general, the re- 
constructed tracks themselves depend on the assumed geometry parameters. This section describes 
how this minimization problem is treated by Millepede [5, 24] — ^the main algorithm used for ITS 
alignment — and presents the first alignment results obtained with cosmic-ray data. 

6.1 General principles of the Millepede algorithm 

Millepede belongs to the global least-squares minimization type of algorithms, which aim at de- 
termining simultaneously all the parameters that minimize the global in Eq. (|6.l|). It assumes 
that, for each of the local coordinates, the residual of a given track f to a specific measured point 
p can be represented in a linearized form as 5; ^ = a - d5t^p/da-\- dt ■ d5t^p/ddt, where a is the 
set of global parameters describing the alignment of the detector (three translations and three ro- 
tations per module) and dt is the set of local parameters of the track. The coiTesponding ;tgiobai 
equation for n tracks with v local parameters per track and for m modules with 6 global parameters 
(N = 6m total global parameters) leads to a huge setof N + vn normal equations. The idea behind 
the Millepede method is to consider the local d parameters as nuisance parameters that are elimi- 
nated using the Banachiewicz identity [25] for partitioned matrices. This allows to build explicitly 
only the set of A'^ normal equations for the global parameters. If needed, linear constraints on the 
global parameters can be added using the Lagrange multipliers. Historically, two versions. Mille- 
pede and Millepede II, were released. The first one was performing the calculation of the residuals, 
the derivatives and the final matrix elements as well as the extraction of the exact solution in one 
single step, keeping all necessary information in computer memory. The large memory and CPU 
time needed to extract the exact solution of a N xN matrix equation effectively limited its use to 
N < 10,000 global (alignment) parameters. This limitation was removed in the second version. 
Millepede II, which builds the matrices (optionally) in sparse format, to save memory space, and 
solves them using advanced iterative methods, much faster than the exact methods. 

6.2 Millepede for the ALICE ITS 

Following the development of Millepede, ALICE had its own implementation of both versions, 
hereafter indicated as MP and MPII, within the AliRoot framework [6]. Both consist of a detector 
independent solver class, responsible for building and solving the matrix equations, and a class 
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Figure 9. (colour online) Left: example of Millepede residuals in the local reference frame of the SPD 
modules before and after the alignment. Right: the corrections to the (p\oc angle obtained in the hierarchical 
SPD alignment with Millepede. 



interfacing the former to specific detectors. While MP closely follows the original algorithm [5], 
MPII has a number of extensions. In addition to the MinRes matrix equation solution algorithm 
offered by the original Millepede II, the more general FGMRES [26] method was added, as well as 
the powerful ILU(k) matrix preconditioners [27]. All the results shown in this work are obtained 
with MPII. 

The track-to-point residuals, used to construct the global x^, are calculated using a parametric 
straight line r{t) = a + bt or helix r{t) = {ox + rcos{t + (po),ay + rsm{t + (po),a^ + t>zt} track 
model, depending on the presence of the magnetic field. The full error matrix of the measured 
points is accounted for in the track fit, while multiple scattering is ignored, since it has no systematic 
effect on the residuals. 

Special attention was paid to the possibility to account for the complex hierarchy of the 
alignable volumes of the ITS, in general leading to better description of the material budget distri- 
bution after alignment. This is achieved by defining explicit parent-daughter relationships between 
the volumes corresponding to mechanical degrees of freedom in the ITS. The aligimient is per- 
formed simultaneously for the volumes on all levels of the hierarchy, e.g. for the SPD the correc- 
tions are obtained in a single step for the sectors, the half-staves within the sectors and the modules 
within the half-staves. Obviously, this leads to a degeneracy of the possible solutions, which should 
be removed by an appropriate set of constraints. We implemented the possibility to constrain either 
the mean or the median of the corrections for the daughter volumes of any parent volume. While 
the former can be applied via Lagrange multipliers directly at the minimization stage, the latter, 
being non-analytical, is applied after the Millepede minimization in a special post-processing step. 
The relative movement 5 of volumes for which the survey data is available (e.g. SDD and SSD 
modules) can be restricted to be within the declared survey precision asurvey by adding a set of 
gaussian constraints 5^/a^^^„^y to the global X'^- 

We report here two example figures to illustrate the bare output results from Millepede II for 
the alignment of the SPD detector (in this case), while the analysis of the alignment quality will be 
presented in the next section. The left-hand panel of Fig. ^ shows an example of the residuals in 
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Figure 10. Left: distribution of Axy|v=o for SPD only, before and after alignment. Right: distribution Axy|v=o 
for track segments reconstructed in the upper and lower parts of SPDh-SSD layers; each track segment is 
required to have four assigned points; SSD survey and Millepede alignment corrections are applied. In both 
cases, the distributions are produced from the sample of events used to obtain the ahgnment corrections. 



SPD (in the local reference frame of the modules) before and after alignment. The right-hand panel 
of Fig. shows the obtained corrections for the (pioc angle (rotation of the volume with respect to 
its zioc axis), indicating that the lai^gest misalignments are at the level of the half-staves with respect 
to the carbon fiber support sectors. 

6.3 Results on alignment quality 

The SPD detector was first aligned using 5 x 10^ cosmic -ray tracks, with two points in the inner 
layer and two points in the outer layer, collected in 2008 with the magnetic field switched off. 
As described in the previous section, the hierarchical alignment procedure consisted in: aligning 
the ten sectors with respect to each other, the twelve half-staves of each sector with respect to the 
sector, and the two modules of each half-stave with respect to the half-stave. 

The following two observables are mainly used to check the quality of the obtained alignment: 
the top half-track to bottom half-track matching at the plane y = 0, and the track-to-point distance 
for the "extra" points in the acceptance overlaps. 

For the first observable, the cosmic-ray track is split into the two track segments that cross 
the upper (y > 0) and lower (y <0) halves of the ITS barrel, and the parameters of the two seg- 
ments are compai^ed at j = 0. The main variable is Ajc3'|,,=o, the track-to-track distance at 3^ = in 
the {x,y) plane transverse to beam line. This observable, that is accessible only with cosmic-ray 
tracks, provides a direct measurement of the resolution on the track transverse impact parameter 
do; namely: CJ/^,|^„(/7t) = V^fydoiPt)- Since the data used for the cun^ent analysis were collected 
without magnetic field, they do not allow us to directly assess the do resolution (this will be the 
subject of a future work). However, also without a momentum measurement, Ajcj| v=o is a powerful 
indicator of the alignment quality, as we show in the following. 

Figure [1^ (left) shows the distribution of A;cy|j,=o for SPD, without and with the alignment 
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corrections. The two track segments are required to have a point in each of the SPD layers and 
to pass, in the transverse plane, within 1 cm from the origin (this cut selects tracks with a similar 
topology to those produced in collisions and rejects tracks that have small incidence angles on 
the inner layer modules). A gaussian fit to the distribution in the range [—100 /xm, + 100 /im] 
gives a centroid compatible with zero and a spread a ^ 50 /^m. For comparison, a spread of 
38 nm is obtained from a Monte Carlo simulation, with the ideal geometry of the ITS (without 
misalignment), of cosmic muons generated according to the momentum spectrum measured by 
the ALICE TPC in cosmic runs with magnetic field. When only the SPD detector is used and the 
tracks are straight lines (no magnetic field), the spread of the Ajc3'| j,=o distribution can be related in a 
simple way to the effective spatial resolution aspatiah which includes the intrinsic sensor resolution 
and of the residual misalignment. For tracks passing close to the beam line (as in our case, with the 
cut at 1 cm), we have: 

2 ^ ^ (''IpDl'^spatiaLSPDl +''SPD2*^spatial,SPD2) ^ « ''sPD 1 + '"SPD2 2 

'^Arvl.-o ~ 7^ ^ ^2 ~ 7Z Z ^2 *^spatial ) \^-^) 

(''SPDl -'^SPD2j (,?^SPD1 -?^SPD2j 

where the inner and outer SPD layers are indicated as SPDl and SPD2, respectively. This relation 
neglects the effect of multiple scattering in the pixels and in the beam pipe, which is certainly one 
of the reasons why the Axyl^^o distribution is not gaussian outside the central region, most likely 
populated by the high-momentum component of the cosmic muons. Using the fit result, CJ^i^,^^ ^ 
50 nm, obtained in the central region [—100 /im, + 100 /^m], we estimate the value aspatiai ~ 14 /^m, 
not far from the intrinsic resolution of about 1 1 /xm extracted from the simulation. However, a 
precise estimation of the effective spatial resolution with this method requires the measurement of 
the track momentum, to account properly for the multiple scattering contribution. The statistics 
collected in 2008 with magnetic field did not allow a momentum-differential analysis. 

The next step in the alignment procedure is the inclusion of the SSD detector As shown 
in section ^ the survey measurements already provide a very precise alignment, with residual 
misalignment levels of less than 5 /^m for modules on the ladder and of about 20 /xm for ladders. 
Because of the limited available statistics (ss 2 x 10^ tracks with four points in SPD and four 
points in SSD), the expected level of alignment obtained with Millepede on single SSD modules is 
significantly worse than the level reached with the survey measurements. For this reason. Millepede 
was used only to align the whole SPD barrel with respect to the SSD baiTcl and to optimize the 
positioning of large sets of SSD modules, namely the upper and lower halves of layers 5 and 6. 
For this last step, the improvement on the global positioning of the SSD layers was verified by 
comparing the position and direction of the pairs of SSD-only track segments built using: two 
points in the upper and two in the lower half-barrel (upper-lower configuration) or two points in 
the inner and two in the outer layer (inner-outer configuration). Before the alignment (only the 
survey corrections applied), the mean of Axy|y=o is (120 ±7) /im and (—1.8 ±0.6) ^um for the 
upper-lower and inner-outer configurations, respectively. After the alignment, it is (—5 ±6) /im 
and (0.5 ±0.6) /^m, respectively, that is, compatible with zero for both configurations. 

The right-hand panel of Fig. [1^ shows the distribution of ATy|y=o for pairs of track segments, 
each reconstructed with two points in SPD and two in SSD, i.e. the merged cosmic-ray track has 
eight points in SPD-i-SSD. It can be seen that, when the SSD survey and the Millepede alignment 
are applied, the distribution is centred at zero and very narrow (FWHM ^ 60 /im), but it shows 
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Figure 11. (colour online) SPD double points in acceptance overlaps. Left: track-to-point Axioc for "extra" 
points before and after alignment. Right: a of the Axioc distributions as a function of the track-to-module 
incidence angle selection; 2008 cosmic-ray data are compared to the simulation with different levels of 
residual misalignment. See text for details. 



non-gaussian tails, most likely due to multiple scattering. A more precise aligmnent of the SSD 
using high-momentum tracks will be performed with the 2009 cosmic-ray and proton-proton data. 

The second alignment quality observable is the Axioc distance between points in the region 
where there is an acceptance overlap between two modules of the same layer. Because of the short 
radial distance between the two overlapping modules (a few mm), the effect of multiple scattering is 
negligible. However, in order to relate the spread of Axioc to the effective resolution, the dependence 
of the intrinsic sensor resolution on the track-to-module incidence angle has to be accounted for. In 
particular, for SPD, due to the geometrical layout of the detector (Fig. ^ left), the track-to-module 
incidence angles in the transverse plane are in general not equal to 90° and they are very different 
for two adjacent overlapping modules crossed by the same track. If Ajcioc is defined as described 
in section ^ the en^or on Ajcioc can be related to the effective spatial resolution of the two points, 

Cspatial) 

f^ivta, = (^^patial ("2) + a,2patial («1 ) COS^ (912) (6.3) 

where the 1 and 2 subscripts indicate the two overlapping points, a, is the incidence angle of the 
track on the module plane, and 912 is the relative angle between the two module planes, which is 
18° and 9° on the inner and outer SPD layer, respectively. Note that, for SSD overlaps on the same 
ladder, we have oti = a2 — 90° and 912 = 0; therefore, cr^noc = \/2o'spatiai> which is the relation we 
used in section 0. 

We start by showing, in Fig. [TT] (left), the track-to-point distance Axioc for the SPD "extra" 
points in the transverse plane, before and after the Millepede ahgnment. The extra points are not 
used in the alignment procedure. The spread of the distribution is a 18 /^m, to be compared to 
a w 15 nm from a Monte Carlo simulation with ideal geometry. An analysis of the Axioc distance 
as a function of the a incidence angle has been performed: five windows on the sum (ai + a2) 
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Figure 12. (colour online) Alignment stability tests. Left: for SPD only, distribution of the A;icy|v=o distance 
obtained when aligning with every second track and checking the alignment with the other tracks. Right: 
track-to-point residuals in the inner SPD layer (track fit in outer SPD layer and in the two SSD layers), for 
B = and B — ±0.5 T (the three histograms are normalized to the same integral). 



of the incidence angles on the two overlapping modules have been considered. These windows 
define increasing ranges of incidence angles from 0° to 50°. Figure [TT] (right) shows the spread 
of the Axioc distribution for the different incidence angle selections: a clear dependence of the 
spread (hence of the spatial resolution) on the incidence angle can be seen. This dependence was 
already observed in SPD test beam measurements [28,29], which were used to tune the detector 
response simulation in the AliRoot software. In the same figure, Monte Carlo simulation results 
are reported for comparison: simulation with ideal geometry (open circles) and with a misaligned 
geometry obtained using a random gaussian residual misalignment (dashed lines: misalignments 
with a = 7 iJ.m and three different seeds; dotted line: misalignments with a = 10 /im). The 2008 
data are well described by the simulation with a random residual misalignment with a « 7 /xm. 
However, this conclusion is based on the assumption that the intrinsic resolution is the same in 
the real detector and in the simulation. Since the intrinsic resolution can slightly vary depending 
on the working conditions of the detector (e.g. the settings used for the bias voltage and for the 
threshold), the value of 7 /xm for the residual misalignment should be taken only as an indication. 
Furthermore, this is an equivalent random misalignment, while the real misalignments are likely 
non-gaussian and to some extent correlated among different modules. 

The robustness of the obtained results was tested by dividing the data sample in two parts 
and using every second track to align the SPD and the others to check the alignment quality. The 
corresponding Ary|y=o distribution is presented in the left-hand panel of Fig. Tl: the distribution is 
centred at zero and has the same a « 50 /xm as in the case of aligning with all tracks. 

Finally, the data with the 0.5 T magnetic field switched on (a few thousand events collected 
at the end of the 2008 cosmic run) were used to perform dedicated checks to evaluate a possible 
effect of the field on the alignment. The alignment corrections extracted from data with B = were 
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applied to data with B = ±0.5 T and alignment quality was verified using both the extra points 
method and the track-to-point residuals method. The distributions of the track-to-point distance for 
extra points in SPD acceptance overlaps for B = 0, +0.5 T and —0.5 T data were found to have 
compatible widths (a [ixm]: 18.3±0.5, 17.8±2.3, 18.4 ± 1.8, respectively) [21]. Another check 
was performed using the track-to-point residuals, calculated by fitting the tracks in the SSD layers 
and the outer SPD layer and evaluating the residuals in the inner SPD layer. In Fig. |T^ (right) the 
comparison between the residuals without magnetic field, with +0.5 T, and with —0.5 T is shown. 
Also in this case, the distributions without field and with the two field polarities are compatible. 

6.4 Prospects for inclusion of SDD in the Millepede procedure 

The alignment of the SDD detectors for the xioc coordinate (reconstructed from the drift time) is 
complicated by the interplay between the geometrical misalignment and the calibration of drift ve- 
locity and to (defined in section The to parameter accounts for the delays between the time 
when a particle crosses the detector and the time when the front-end chips receive the trigger signal. 
Two methods have been developed in order to obtain a first estimate of the to parameter. The first, 
and simpler, method consists in extracting the to from the minimum measured drift time on a large 
statistics of reconstructed SDD points. The sharp rising part of the distribution of measured drift 
times is fitted with an error function. The to value is then calculated from the fit parameters. The 
second method measures the to from the distributions of residuals along the drift direction (xioc) be- 
tween tracks fitted in SPD and SSD layers and the corresponding points reconstructed in the SDD. 
These distributions, in case of miscalibrated to, show two opposite-signed peaks corresponding to 
the two separated drift regions of each SDD module, where electrons move in opposite directions 
(see Fig. ^ right). The to can be calculated from the distance of the two peaks and the drift ve- 
locity. This second procedure has the advantage of requiring smaller statistics, because it profits 
from all the reconstructed tracks, with the drawback of relying on SDD calibration parameters (the 
drift velocity and possibly the correction maps). Moreover, being based on track reconstruction, it 
might be biased by SPD and/or SSD misalignments. 

Depending on the available statistics, the to determination with these two methods can be done 
at the level of SDD ban^el, SDD ladders or SDD modules. The to parameter needs actually to be cal- 
ibrated individually for each of the 260 SDD modules, because of differences in the overall length 
of the cables connecting the DAQ cards and the front-end electronics. In particular, a significant 
difference is expected between modules of the A (z > 0) and C sides (z < 0), due to the « 6 m 
difference in the length of the optical fibres connecting the ITS ladders to the DAQ cards. With the 
first 2000 tracks, it is possible to determine the to from track-to-point residuals for 4 sub-samples 
of modules, i.e. sepai^ating sensors connected to sides A and C of layers 3 and 4. An example of 
residual distributions for the left and right drift sides of the modules of layer 4 side C is shown 
in Fig. |T3[ The Millepede alignment corrections for SPD and SSD are applied in this case, and it 
has been checked that, if they are not applied, the centroid positions in this figure are not affected 
significantly, while the spread of the distributions increases, as it could be expected. A difference 
of 25 ns between sides A and C of each SDD layer has been observed, in agreement with the 6 m 
difference in fibre lengths (the propagation time of light in optical fibres is 4.89 ns/m). With larger 
statistics (35,000 tracks), it is possible to extract the to for each half-ladder, which requires pro- 
ducing 36(ladders)x2(A/C sides) pairs of histograms like the ones shown in Fig. [T3| (left). Further 
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track-to-point residuals in x^^^ [cm] drift coordinate x,^^ [cm] 

Figure 13. SDD calibration and alignment. Left: distribution of track-to-point residuals in the two drift 
regions for the SDD modules of layer 4 side C (z < 0); tracks are fitted using only their associated points 
in SPD and SSD; the Millepede alignment corrections for SPD and SSD are included, as well as the SSD 
survey. Right: residuals along the drift coordinate for one SDD module as a function of drift coordinate after 
Millepede alignment with only geometrical parameters and with geometrical-ncalibration parameters. 



cable-length differences, which introduce to difference, exist at the level of individual half ladders 
and at the level of individual modules. These differences are of the order of 1.5 m (7 ns) and 20 cm 
(< 1 ns) respectively, but have not been measured in detail yet, because of the limited size and 
coverage of the cosmic -ray tracks data sample. It should be noted that given the 6.5 /im/ns of 
drift velocity, a bias of 1 ns on the to can lead to a significant effect on the reconstructed position 
along the drift coordinate xioc. 

After a first calibration with these methods, a refinement of the to determination is obtained by 
running the Millepede minimization with the to as a free global parameter for each of the 260 SDD 
modules. Similarly, the drift velocity is considered as a free parameter for those SDD modules 
(about 35%) with mal-functioning injectors. For these modules a single value of drift velocity 
is extracted for the full data sample analyzed, thus neglecting the possible dependence of drift 
velocity on time due to temperature instabilities. However, this is not a major concern since the 
drift velocity was observed to be remarkably stable during the data taking [17]. This allows to 
assess at the same time geometrical alignment and calibration parameters of the SDD detectors. 
About 500 tracks are required to align and calibrate a single SDD module. An example is shown 
for a specific SDD module in the right-hand panel of Fig. |T3|, where the xioc residuals along the 
drift direction are shown as a function of xioc- The result obtained using only the geometrical 
rotations and translations as free parameters in the Millepede minimization is shown by the circle 
markers. The clear systematic shift between the two drift regions (xioc < and .X]oc > 0) is due to 
both miscalibrated to and biased drift velocity (this is a module with non-working injectors). These 
systematic effects are no longer present when also the calibration parameters are fitted by Millepede 
(square markers). It should be pointed out that the width of the SDD residual distributions shown 
in Fig. [T3| does not correspond to the expected resolution on SDD points along drift coordinate 
because of jitter between the time when the muon crosses the detectors and the SPD FastOR trigger, 
which has an integration time of 100 ns. For the about 100 SDD modules with highest occupancy. 
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the statistics collected in the 2008 cosmic run allowed to check the reliability of the calibration 
parameters {to and drift velocity) extracted with Millepede by comparing the values obtained from 
independent analyses of two sub-samples of tracks. From this study, a precision of 0.025 /^m/ns 
for the drift velocity and 10 ns for to was estimated. It should be noted that these precisions are 
limited by the available statistics as well as by the trigger jitter effect mentioned above. 



7. SPD alignment with an iterative local method 

We developed an alignment method that performs a (local) minimization for each single mod- 
ule and accounts for correlations between modules by iterating the procedure until convergence 
is reached. A similar approach is considered by both the CMS and ATLAS experiments [30-32]. 
The main difference between this method and the Millepede algorithm is that only in the latter 
the correlations between the alignment parameters of all modules are explicitly taken into account. 
Conversely, the local module-by-module algorithm assumes that the misalignments of the modules 
crossed by a given track are uncorrected and performs the minimization of the residuals indepen- 
dently for each module. The comparison of the alignment parameters from this method and from 
Millepede would provide a further validation of the results achieved with the Millepede. 

In the local method we minimize, module-by-module, the following local function of the 
alignment parameters of a single module: 



^local(^tra,Arot) — ^ ^t^p^ t ,p ,p 
tracks 

52 (j, - ArotTp - atra)^ (Vj + Vp)"' (r, - Arot^p - fltra) • (7. 1) 

tracks 



Here, the sum runs over the tracks passing through the module, fp is the position of the measured 
point on the module while is the crossing point on the module plane of the track t fitted with 
all points but Tp. \t and Vp are the covariance matrices of the crossing point and of the measured 
point, respectively. The six alignment parameters enter this formula in the vector Otra; the alignment 
correction for the position of the centre of the module, and in the rotation matrix Ajot, the alignment 
correction for the orientation of the plane of the module. The alignment correction is supposed to 
be small so that the rotation matrix can be approximated as the unity matrix plus a matrix linear 
in the angles. In this way, the is a quadratic expression of the alignment parameters and the 



minimization can be performed by simple inversion. The X\oca\ function in Eq. (7. 1 ) can be written 
in the same way also for a set of modules considered as a rigid block. The track parameters are not 
affected by the misalignment of the module under study, because the track point on this module 
is not used in the fit, while the positions of the crossing points are affected, because the tracks are 
propagated to the plane of the module defined in the ideal geometry. This is taken into account by 
adding a large error along the track direction to the covariance matrix of the crossing point. 

Given that this is a local method, it is expected to work best if two conditions are fulfilled: 
the correlation between the misalignments of different modules is small and the tracks used to 
align a given module cross several other modules. In order to limit the bias that can be introduced 
by modules with low statistics, for which the second condition is normally not met, we align the 
modules following a sequence of decreasing number of points. To reduce the residual correlation 



-30- 




Figure 14. SPD alignment quality results for the iterative local method. Left: track-to-track Axy|y=o distri- 
bution defined in section |6.3| , using only SPD points. Right: track-to-point Axioc distribution for extra points 
in acceptance overlaps. 



between the alignment parameters obtained for the different modules, we iterate the procedure 
until the parameters converge. Simulation studies with misalignments of the order of 100 /xm have 
shown that the convergence is reached after about 10 iterations. 

For the ITS alignment using the 2008 cosmic-ray data, we aligned only the SPD modules 
using this method. Like for Millepede, we adopted a hierarchical approach. Given the excellent 
precision of the SSD survey measurements, we used these two layers as a reference. We aligned 
as a first step the whole SPD barrel with respect to the SSD, then the two half-barrels with respect 
to the SSD, then the SPD sectors with respect to the SSD. In the last step, we used SPD and 
SSD points to fit the tracks and we aligned the individual sensor modules of the SPD. Figure 
shows the top-bottom track-to-track Ajcj^jy^o distribution obtained using only the SPD points (left- 
hand panel) and the track-to-point Axioc for the double points in acceptance overlaps (right-hand 
panel), after alignment. Both distributions are compatible (mean and sigma from a gaussian fit) 
with the corresponding distributions after Millepede alignment. This is an important independent 
verification of the Millepede results. Since the two methods are in many aspects independent, 
comparing the two sets of alignment parameters could provide a check for the presence of possible 
systematic trends. Figure [T3| shows the correlation of the inner SPD layer parameter values obtained 
with the iterative method and those obtained with Millepede. A correction was applied to account 
for a possible global roto-translation of the whole ITS, which does not affect the quality of the 
alignment and can be different for the two methods. The closed (open) markers represent the 
modules with more (less) than 500 track points. Most of the modules are clustered along the 
diagonal lines where the parameters from the two methods are exactly the same. There are some 
outlier modules, that are far from the ideal result. However, these outliers mostly correspond to 
modules with low statistics (open markers) at the sides of the SPD barrel. 
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Figure 15. Correlation between the alignment parameters obtained from Millepede (horizontal axis) and the 
iterative method (vertical axis), for the inner SPD layer modules. Modules with more (less) than 500 points 
are represented by the closed (open) markers. 



8. Conclusions 

The results on the first alignment of the ALICE Inner Tracking System with cosmic-ray tracks, 
collected in 2008 in the absence of magnetic field, have been presented. 

The initial step of the alignment procedure consisted of the validation of the survey measure- 
ments for the Silicon Strip Detector (SSD). The three methods applied for this purpose indicate 
that the residual misalignment spread for modules on ladders is less than 5 nm, i.e. negligible with 
respect to the intrinsic resolution of this detector in the most precise direction, while the residual 
misalignment spread for the ladders with respect to the support cones amounts to about 15 jim. 

The procedure continues with track-based software alignment performing residuals minimiza- 
tion. We presented the results obtained with a sample of about 10^ cosmic-ray tracks, reconstructed 
in events selected by the FastOR trigger of the Silicon Pixel Detector (SPD). We mainly use the 
Millepede algorithm, which minimizes a global of residuals for all alignable volumes and a 
large set of tracks. 

We start from the SPD, which is aligned in a hierarchical approach, from the largest mechanical 
structures (10 support sectors) to the 240 single sensor modules. About 90% of the latter were active 
during the 2008 cosmic run, and about 85% had enough space points (> 50) to perform alignment. 
Then, we align the SPD barrel with respect to the SSD barrel. The SSD coverage provided by the 
cosmic -ray tracks is insufficient to align the SSD at the level of ladders, especially for the ladders 
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close to the horizontal plane y = 0. Therefore, for the time being we only align the SSD at the level 

of large sets of ladders. 

The two intermediate ITS layers, the Silicon Drift Detectors (SDD), represent a special case, 
because the reconstruction of one of the two local coordinates requires dedicated calibration pro- 
cedures (drift velocity and drift time zero extraction), which are closely related to the aUgnment. 
Indeed, one of the approaches that we are developing for the time zero calibration is based on the 
analysis of track residuals in a standalone procedure, initially, and then directly within the Mille- 
pede algorithm. Once these procedures are stable and robust, the SDD will be included in the 
standard alignment chain. For all six layers, the completion of the alignment for all modules will 
require tracks from proton-proton collisions; a few 10^ events (collected in a few days) should 
allow us to reach a uniform aUgnment level, close to the target, over the entire detector. 

We use mainly two observables to assess the quality of the obtained alignment: the matching 
of the two half-tracks produced by a cosmic-ray particle in the upper and lower halves of the ITS 
barrel, and the residuals between double points produced in the geometrical overlaps between ad- 
jacent modules. For the SPD, both observables indicate an effective space point resolution of about 
14 [xm in the most precise direction, only 25% worse than the resolution of about 1 1 [xm extracted 
from the Monte Carlo simulation without misaUgnments. In addition, the measured incidence angle 
dependence of the spread of the double points residuals is well reproduced by Monte Carlo simu- 
lations that include random residual misalignments with a gaussian sigma of about 7 /xm. Further 
confidence on the robustness of the results is provided, to some extent, by the cross-checks we per- 
formed using a small data set with magnetic field switched on and, mainly, by the comparison of the 
Millepede results to those from a second, independent, alignment method. This second method, 
which iteratively minimizes a set of local module-by-module functions, yields, compared to 
Millepede, a similar alignment quaUty and a quite compatible set of alignment corrections. 

Using the present data set with magnetic field off, since the track momenta are not known, 
the multiple scattering effect, which is certainly not negligible, cannot be disentangled from the 
residual misalignment effect. Therefore, a more conclusive statement on the SPD residual mis- 
alignment will be possible only after the analysis of cosmic-ray data collected with magnetic field 
switched on. The same apphes for combined tracking with SPD, SDD and SSD points: in this case, 
the momentum-differential analysis of the transverse distance between the two half-tracks (upper 
and lower half-barrels) will allow us to measure the track transverse impact parameter resolution, 
which is a key performance figure in view of the ALICE heavy flavour physics program. 
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